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Abstract 

There are problems in physics and particularly in field theory which 
are defined by complex valued weight functions e~ s where S is a poly- 
nomial action S : IR n — > C. The conditions under which a convergent 
complex Langevin calculation correctly simulates such integrals are 
discussed. All conditions on the process which are used to prove proper 
convergence are defined in the stationary limit. 
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1 Introduction 



Complex Langevin (CL) methods have turned out to be quite useful in the 
calculation (simulation) of high dimensional integrals over complex valued 
weight functions of the form e~ s , where S is the action or the Hamiltonian 
of some physical system. Since there is no formal restriction to a real valued 
drift term for Langevin equations, the application of CL is convincingly sim- 
ple 0. Unfortunately one has to deal with two problems of uncertainty. The 
first is that it is apriori unknown whether the process will converge at all. 
The second problem is that, although the process has converged, it will not 
necessarily give the correct answer. This is, that long time averages of such 
a process do not necessarily simulate the complex valued weight function 
integrals. CL is known to sometimes give the wrong answer (see e.g. ||). 

Several attempts have been made to understand CL (e.g. see references 
H, H|, |J). For some simple actions the behavior of CL can be improved by 
modifying the drift term with an appropriate kernel, but for general problems 
the choice of the kernel is not clear ||. Recently progress has been made in 
the comprehension of the results which one gets from a convergent process 
[§, [7|. In particular the assumptions needed to guarantee correct results for 
convergent processes on certain compact manifolds (Si, Sz) turn out to be 
surprisingly simple and easy to verify in a numerical simulation. Contrary 
to that, many assumptions are used to prove the behavior of processes on 
M n driven by polynomial actions and moreover these assumptions are rather 
technical ||. 

For polynomial actions a lot of attention has been given to the existence 
of a pseudo Fokker-Planck (F-P) equation which describes the dynamics of 
a possibly equivalent complex valued weight function p|. In earlier inves- 
tigations especially the spectrum of this operator played a major role [p]]. 
But statements on the properties of the spectrum are not sufficient to draw 
conclusions on the correctness or the convergence of CL ||. Certainly, if one 
can show that the pseudo F-P equation exists and that the real part of the 
spectrum of the operator is semidefinite then CL converges but not necessar- 
ily to the desired result. Further conditions must hold (see @). Except for 
very simple cases it is hard and most unlikely to get exact information on the 
complete spectrum. Certainly there always exist the real F-P operator for 
the process and the convergence of the process follows if one can prove that 
the operator has a unique nonnegative integrable solution to the zero eigen- 
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value. But this is also very hard and so far there is no classification scheme 
for actions which have the suitable properties. So, to get information on the 
convergence for any problem one must check either the existence and the 
whole spectrum of the pseudo F-P operator or the zero eigenvalue properties 
of the real F-P operator. In practice therefore the question of convergence 
still remains a matter of experiment and experience. 



2 Proper Convergence 

Let us now turn to the main purpose of the paper and examine in a rigorous 
fashion the conditions under which CL, if convergent, gives the right answer. 
To demonstrate this, several conditions at finite time have been put on the 
process in reference |§. In this approach fewer conditions on the process 
are used and these conditions are put forward to t —>■ oo. For simplicity the 
discussion and the formulas are restricted to the one dimensional case. All 
following statements allow for an immediate generalization to lR n . It will be 
assumed that the system of interest is described by a complex polynomial 
action of degree iV 

N 

S(x) = £a n x" . (2.1) 

n=0 

S : R -> C such that e~ s G S(M). S(R) is the Schwartz space of C°° 
functions of rapid decrease. With g(x) a polynomial of degree M it is thus 
guaranteed, that the quantities of physical interest 

(g(x)) = ^ j^g(x)e' s ^dx, (2.2) 
Af = [ e~ s(x) dx , (2.3) 



do exist, provided < \J\f\. If S would be real valued everything would be 
now straight forward ergodic theory and the longtime averages computed 
with the Langevin equation would reproduce the ensemble average of the 
system. ||) 

In the complex case analytic continuation leads to the following stochastic 
differential equation. 

dZ(t) = F(Z(t))dt + dW(t) , (2.4) 
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with the drift term 

W(t) is a standard Wiener process with zero mean and covariance 

E {Wit^Wih)) = mmih, t 2 ). (2.6) 



Equation |274] is the so called CL equation. This equation has a locally unique 



solution which is defined up to a random explosion time M. In particular 



equation |2.4j describes a two dimensional diffusion process. 

dX(r) = G(X(t), Y(r))dr + dW(r) , (2.7) 

dY(r) = H(X(r),Y(T))dT . (2.8) 
With S(z) = u(x, y) + iv(x, y), we have 

G{X , V) = -\ ! ^A > H(x , y) = l^lA. ( 2 . 9 ) 



Special for this process is that equation ETBI looks like a deterministic equation 



due to the zero diffusion coefficient. Nevertheless this is a stochastic equation 
through the dependence on X(t). The singular diffusion matrix causes a lot 
of problems. So, contrary to the real action case it is in general not possible 
to determine from the drift and diffusion terms whether there exists a unique 
stationary distribution density for this process ||. As already mentioned in 
the introduction there is no general proof on the existence of a stationary 
distribution density. For the moment let us assume that for the process 
X(t),Y(t) there exists a unique stationary distribution density f(x,y). The 
idea behind CL then is that 



Jim E(g(X(t) + iY(t))) 



t— >oo 



g(x + iy)f{x, y)dxdy = jj= L 9( x ) e S(x) dx. 



(2.10) 



might hold. 
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Assume: 



1. S is a complex valued polynomial action of degree N such that 



and 



2. For 



-S( x ) d 



x 



> 0. 



j(M) = E{e ikZ{t) ) = J^ 2 e lk{x+iy) f{x,y,t)dxdy 



the limit t —>■ oo exists pointwise and 

lim c T (k) = Cqo (A;) £ «S( 



3. Further 



lim 

t— >oo 



£(Z n (t)e 



ikZit)- 



< oo for all < n < N - 1, A; G 



(2.11) 
(2.12) 

(2.13) 
(2.14) 

(2.15) 



Equation |2.10| then holds at least for g(z) a polynomial of degree M < N — 1. 
Moreover equation p.!0| holds for any higher moment E(Z n (t)), n > N which 
exist for t — > oo. 

From assumption 2 we know that there is a to < oo such that c(k, t) exists 
and from assumption 3 that there is a ti < oo such that E{Z n {t)e ikm ) exists. 
Applying the Ito rule one gets with F(z) as defined in pT5| 



a^(e 



ikE (e ikz ®F(Z(t))) - —E ( 



dt 



,ikZ(t) 



(2.16) 



Due to assumptions 2 and 3 equation [2.16| exists for t' = max(t ,ti)- As a 
side result we get that, if c(k,t) G C 7V_1 (M) with repsect to k, equation |2.16 
can be understood as the dynamical equation for c(k,t). 



dc(k,t) 
dt 



2 n t! 1 dk 



(~ti r ) n - 1 c(k,t)-—c(k,t) 



(2.17) 



Note that if assumption 3 does not hold, equation |2.16| can also not be defined 
in the sense of distributions. This is because we are not simply dealing with 
Fourier transforms but with their possibly not existing analytic continuations. 
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Let us define now h(x) as 
h(x) = 



1 

2^ 



Coo {k)e- lkx dk . 



(2.18) 



From assumption 2 follows that h(x) G <S(K). Using equation |2.18| and 
assumption 3 

(2.19) 



lim E(Z n (t)e 



t— >oo 



x n e h(x)dx 



for < n < N — 1 and fcel. Applying the above result to |2.16| one obtains 
in the limit t — > oo 



Q = ik e tkx F(x)h(x)dx - 



e ikx h{x)dx 



(2.20) 



Integrating the right hand side of equation [2.20| by parts gives that h(x) is a 
L 1 (R, dx) zero eigenvalue solution of a F-P type differential operator with a 
complex drift term (pseudo F-P operator). 



ld_ 

2dx 



dS(x) | d 
dx dx 



h(x) = Th(x) = 0. 



T has two zero eigenvalue solutions. One is 

h t {x) ~ e~ s{x) G S(R) 



(2.21) 



(2.22) 



which fits to assumption 2, since as the Fourier transform of c^k) it must 
be a Schwartz function. For the second solution 



one can show that 



h 2 (x) 



h 2 {x) = 0{- 



-S(x) 



' {y) dl 



■''() 



X 



N-l 



for \x\ — > oo 



(2.23) 



(2.24) 



This contradicts assumption 2. So, the only possible solution is the one 
proportional to e and 



1 



lim E(Z n (t)e lkZ{t) ) - , 



x n e ikx e s ^dx 



(2.25) 
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for < n < N - 1 and k G R. If further E(Z n (t)), n > N for t -> oo exist 
then 

rf" Coo (A;) 



lim E(Z n (t)) 



dk T 



f x n e S{x) dx _ 26 N 



Let us now briefly discuss the assumptions. Polynomial actions are very 
natural since most physical systems defined on W 1 have polynomial actions. 
Since these actions must be bounded from below it follows that e~ s G S. 
Condition 2 must be there otherwise the solution h 2 (x) cannot be excluded. 
With the correctness requirement on CL that 

lim E(e ikm ) = \- r \ e lkx e- s{x) dx (2.27) 

t^oo v ' Mm 

this condition is also a necessary condition. Assumption 3 looks technical, 
but is so far required to relate hi(x), the Fourier transform of c^k), to the 
Fokker-Planck type operator T. This condition finally allows to show the 
correctness of CL. It would be nice to eliminate assumption 3 by showing 
that it follows from assumption 2. Unfortunately the integral transform de- 
fined by equation [2.13| is not an injective mapping. To the authors knowledge 
the nature of this integral transform has not been analized in the literature. 
At present, without more detailed information on the probability density (in 
general not available), it is perhaps impossible to draw conclusion on the 
properties of the function from the properties of its image. In a numeri- 
cal simulation certainly such mathematical criteria a hard to verify exactly. 
Nevertheless experience tells us that when plotting such expectation values 
(coo(fc), E(Z n e lkZ )) one gets a very clear sign of the quality of the result [110 



3 Conclusions 

The criteria under which a convergent CL simulation leads to correct results 
have been significantly simplified. The assumptions used in the present proof 
are much closer to a numerical verification than the one used in reference 
|J. Unfortunately a complete theory of CL is still lacking. However the 
situation that it was generally neither apriori nor aposteriori possible to 
prove convergence to the desired result has been ameliorated in as far as a 
simple aposteriori proof is now possible. 
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